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Glassy dynamics: effective temperatures and intermittencies from a two-state model 
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We show the existence of intermittent dynamics in one of the simplest model of a glassy system: 
the two-state model, which has been used— to explain the origin of the violation of the fluctuation- 
dissipation theorem. The dynamics is analyzed through a Langevin equation for the evolution of the 
state of the system through its energy landscape. The results obtained concerning the violation factor 
and the non-Gaussian nature of the fluctuations are in good qualitative agreement with experiments 
measuring the effective temperature and the voltage fluctuations in gels and in polymer glasses. The 
method proposed can be useful to study the dynamics of other slow relaxation systems in which 
non-Gaussian fluctuations have been observed. 



I. INTRODUCTION 

Complex systems are often distinguished by the exis- 
tence of a very intricate free energy landscape consist- 
ing of many barriers which the system has to overcome 
to evolve. It is precisely the presence of these barriers 
the responsible for the slow relaxation dynamics which 
manifests in the appearance of peculiar phenomena as 
aging, lack of a fluctuation-dissipation theorem and in- 
termittencies caused by the presence of large fluctuations. 
These features, predicted and observed in systems of dif- 
ferent nature as glasses, granular flows, foams, crum- 
pled materials and in the dynamics of the disordered 



system* 
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have attracted the interest of many re- 



searchers during the last years with the purpose of de- 
scribing the main features of slow relaxation dynamics^. 

In a previous papesi, we have proposed a minimal re- 
laxation model aimed at characterizing the dynamics of 
a system relaxing in two very different time scales, which 
are related to inter-well and intra-well relaxation pro- 
cesses. Two scenarios were analyzed for this purpose. 
In the first of them, the system may explore the whole 
reaction coordinate space undergoing a diffusion process 
described by a Fokker-Planck equation^, which accounts 
for the intra-well and inter-well relaxations. In the other, 
obtained from the first one by eliminating the fast vari- 
able, the system undergoes an activated process. In spite 
of its simplicity the model shows some of the peculiar 
features of the dynamics of slow relaxation systems and 
proposes an explanation of why and how the fiuctuation- 
dissipation is violated. It was found that the violation 
factor or effective temperature depends on the observ- 
able and on the initial populations in the wells. This 
result shows that the effective temperature does not uni- 
vocally characterize the thermal state of a glassy system 
undergoing activated dynamicsifi. 

Our purpose in this paper is to use that model to ex- 
plain the presence of intermittencies in the dynamics of 



a system in a glassy phase and the non-Gaussian nature 
of the probability distribution function, which have re- 
cently been observed in measurements of the dielectric 
properties of gels and polymer glasse a--'-'" !, and in some 
theoretical studies of spin-glass modelsi^^i^. 

The paper is organized as follows. In Section 2, we 
analyze some of the main traits of glassy dynamics from 
a two-state model. Section 3 is devoted to present the 
results concerning the intermittent behavior and the non- 
Gaussian nature of the fluctuations. Some conclusions 
and perspectives are presented in the final section. 



II. GLASSY DYNAMICS FROM A TWO-STATE 
MODEL 

In a two-state model, the minimal relaxation model for 
glassy systemaiSiiSiii, one assumes that the process con- 
sists of two main steps: a slow relaxation, in which the 
coordinate characterizing the state of the system jumps 
from a potential well to the next one, and a fast equi- 
libration process in the well. It has been showni that 
the dynamics of the system can be analyzed in terms of 
a Fokker-Planck equation describing a diffusion process 
through the free energy landscape <l>(7)iS, 



dt 
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(1) 
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In the simplest case, $(7) would be just a bistable 
potential. In the previous equation, ^(7, t) is the prob- 
ability distribution function which depends on the order 
parameter or reaction coordinate 7, D is the diffusion 
coefficient, T is the temperature of the bath and /cb the 
Boltzmann constant. When the height of the barrier sep- 
arating the two minima of the potential is large enough 
compared to thermal energy the systems achieves a state 
of quasi-equilibrium in each well. The evolution of the 
system then proceeds by jumps from one well to the other 
undergoing an activated process. The dynamics corre- 
sponding to this situation can be obtained by eliminat- 
ing the fast degrees of freedom in such a way that it can 



2 



be characterized simply by the populations at each well. 
The Fokker-Planck equation then reduces to the follow- 
ing kinetic equations governing the population dynamics 
at both wellsiiia 



dni 
~dt 



dn2 
~dt 



= -j{t)-f{t). 



where the current j{t) is defined as: 



jit) =3^ - JV 



k^7ii — k^n2 



(2) 



(3) 



Here ni{t) and n2{i) are the populations at each well 
and k^^^ are the forward and backward reaction rates 
given by 



D^^"{ni.2W{lo)\ 

2TTkRT 



exp 



$(71,2) - $(70) 
knT 



(4) 

where 71 and 72 denote the position of the two minima 
and 7o the position of the top at the barrier. The noise 
term j'' [t) results from the coarsening of the correspond- 
ing random term J^('^,t) in the Langevin equation re- 
lated to the Fokker-Planck equation QJ^i. Whereas the 
latter satisfies the fiuctuation-dissipation theorem inher- 
ent to the diffusion process along the reaction coordinate 
formulated as^fi 
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FIG. 1: The effective temperature plotted as a function of 
the frequency (in arbitrary units) for different values of the 
waiting time. 



e{t) — EoOit — t^) is a generelized force that is coupled 
to the observable). The response is then given by 



R{t, 



dim)) 



(7) 
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where the average of the observable is defined as 



( J'-(7, t)r{^', t')) = 2Z?(p(7, i))<5(7 - ^')5it - t') (5) 

the former does not obey a similar expression. Its cor- 
relation is given bjii^: 



{J^mit')) - {k^{ni)+k^{n2))5{t-t') ^ 2k^nlH{t~t') 

(6) 

where the last term is the value of the correlation at 
equilibrium, and n!^ is the equilibrium population in the 
first well. 

This result clearly indicates that the violation of the 
fiuctuation- dissipation theorem is precisely due to the 
coarsening of the description. When one eliminates 
the fast variables reducing the dynamics to that of an 
activated process the system cannot progressively pass 
through local equilibrium states from one well to the 
other but proceeds by jumps and consequently remains 
outside equilibrium. This explains why the fiuctuation- 
dissipation theorem, a result strictly valid when the fiuc- 
tuations take place around an equilibrium stateil, is not 
fulfilled. 

From the model proposed we can analyze the non- 
equilibrium response of the system. Let us assume an 
external perturbation —e{t)0{t) that is plugged in at 
the waiting time t^,, defined as the time elapsed af- 
ter quenching {0{t) is the observable considered and 



{50{t)) 



Oij)Sp{j,t)dj^iO,-02)Sni 



(8) 



whereas the values Oi and O2 are the values of the ob- 
servable at the minima 1 and 2, respectively. On the 
other hand, the correlation is 



Co(i,i») =e-(*-*™)/"(Oi-02)[j^(Oi-Oo)-i^(02-Oo)]. 

(9) 

Both quantities satisfy the relation 



Ro{t, tw) 



1 d 



Co{t, tw), 



(10) 



in which the quantity T^^j plays the role of an "effective" 
temperature given by 



(11) 



Here t ^ = {k^ + k^) is the relaxation time of the 
process and the parameter A has the form 



_ fc^(ni(0))(Oi - Oo) - fc^(n2(0)(O2 - Op) ^^^^ 



k^nl\0i^02) 
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Since the equilibrium state is achieved for — > oo, one 
can easily verify that the effective temperature coincides 
with that of the bath and the relation ifTTijl becomes the 
fluctuation-dissipation theorem. 

An important consequence of our previous analysis is 
that the effective temperature is not a robust quantity 
since it depends on the observable considered as well as 
on the initial populations at the wells. We then conclude 
that in the case of an activated process that quantity 
is a parameter measuring the distance of the system to 
the equilibrium state but it is not universal and conse- 
quently can hardly be considered as a thermodynamic 
temperatureii. 

Our expression En. lfTTll can be used to reproduce the 
dependence of the violation factor on the frequency ob- 
served in experimentsiiii^. If we assume that the whole 
relaxation process takes place through consecutive acti- 
vated processes for which Eq.JJJ applies, we can infer a 
global behavior of that quantity by identifying the in- 
verse of the relaxation time with a frequency. By consid- 
ering the case A = 0, in which the violation factor does 
not depend on the observable , we then conclude that 
the effective temperature at low frequencies behaves as: 
Teff ~ Lu~^. In the experiments one flnds that the effec- 
tive temperature decreases when increasing the frequency 
following a power law whose exponent is in between —1.1 
and —1.2,^^. As in the experiments, it is found that the 
effective temperature tends to the bath temperature for 
very large waiting times (see Fig. 1). 



and 72 = 3 and the parameters ai and 02 have been 
chosen arbitrarily as oi = 6 and 02 = 9, ensuring the 
continuity of the observable. The location of the maxi- 
mum of the observable has been placed at 7 = 1, to show 
in a clearer way the occurrence of a transition between 
the two wells. When the state is at the top of the barrier 
the probability to go back to the minimum 71 is 50%. 

The value of a should be consistent with the non- 
stationary nature of the process and has to be intrin- 
sically related to the waiting time. Therefore, we have 
use different values of this parameter to mimic different 
values of the waiting time. 

In Fig. Illll we have represented the value of the observ- 
able 0(7) corresponding to a single trajectory obtained 
from the simulations of Ea. ljl3|l . with a = 1.5. We ob- 
serve the presence of an intermittent event, which is the 
signature of a jump from one well to another. Therefore, 
the presence of intermittencies is a consequence of the 
activated nature of the relaxation process. 

Averaging over many trajectories we have obtained 
the probability distribution function for the observable, 
which has been plotted in Fig. 4, for different values of 
a. The non-Gaussian form of the curves is, in our model, 
a consequence of the non-parabolic form of the poten- 
tial. Intermittencies add more weight to the tails of the 
distribution, thus stressing its non-Gaussian nature. 



IV. CONCLUSIONS 



III. INTERMITTENT DYNAMICS 

We will analyze in this section the intermittent behav- 
ior of the fluctuations taking place when the system is 
quenched below the glass transition. This type of be- 
havior has been reported in experiments measuring the 
fluctuation spectrum and the response of a Laponite solu- 
tion and of a polymer glassii^. To this purpose, we will 
particularize the model discussed in Sec. 2 to the case 
of the quartic potential plotted in Fig. 2. The dynamics 
of this model is governed by the Fokker-Planck equation 
Q or through the equivalent Langevin equation 

J = -7(7-7i)(7-72) + J''. (13) 

where has been considered a Gaussian white noise 
process having, as a flrst approximation, a constant am- 
plitude: D{p{j,t)) « a. We will see that this approxi- 
mation is enough to explain the experimental results. 

By numerical simulation of the Langevin equation II13II . 
we have computed the value of the following observable 

0(7) = I ""1^"^^! (i4) 
^ " [ 02(7 - 72) if 7 > 1 ' ' 

which has a piece- wise dependence on 7, having a zero 
value at 71 and 72. In our case we have taken 71 = —2 



In this paper we have used the activation over a bar- 
rier to model the evolution of a system when is quenched 
below the glass transition. Using as a model a simple 
quartic potential, we have proved the existence of an in- 
termittent dynamics in the aging process caused by the 
presence of large fluctuations. Similar behavior has been 
observed in recent experiments. We have shown that the 
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FIG. 2: Quartic potential used in our model. 
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FIG. 3: The evolution of a single realization of the observable 
0(7), corresponding to q = 1.5. The single step-time has 
been taken as 0.01. We have considered 10000 time intervals. 



probability distribution function corresponding to those 
events deviates from its Gaussian form observed for sys- 
tems close to equilibrium and exhibits exponential tails 
as those encountered experimentally. 

The model proposed and its possible generaliza- 
tions could, with the help of the stochastic processes 
theory22iSi, constitute a useful tool to characterize the 
dynamics of systems far from equilibrium for which a 
common phenomenology in the behavior of the fluctua- 



tions is being found. 




FIG. 4: Probability distribution function for different values 
of Of as a function of the values of the observable, obtained by 
performing simulation of En. 11811 . The solid line is a Gaussian 
curve plotted just for comparison. 
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